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ff^ ' Abstract 

^ ^ Recent progress on QCD sum rule determinations of the light and heavy quark masses is reported. In 

I , the light quark sector a major breakthrough has been made recently in connection with the historical 

systematic uncertainties due to a lack of experimental information on the pseudoscalar resonance spectral 
(— I , functions. It is now possible to suppress this contribution to the 1% level by using suitable integration 

kernels in Finite Energy QCD sum rules. This allows to determine the up-, down-, and strange-quark 
masses with an unprecedented precision of some 8 — 10%. Further reduction of this uncertainty will be 
possible with improved accuracy in the strong coupling, now the main source of error. In the heavy quark 
. sector, the availability of experimental data in the vector channel, and the use of suitable multipurpose 

' integration kernels allows to increase the accuracy of the charm- and bottom-quarks masses to the 1% 

00 ; level. 
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1 Introduction 



Quark and gluon confinement in Quantum Chromodynamics (QCD) precludes direct experimental mea- 
surements of the fundamental QCD parameters, i.e. the strong interaction coupling and the quark masses. 



- 1 — I 

X 

, Hence, in order to determine these parameters analytically one needs to relate them to experimentally 

measurable quantities. Alternatively, simulations of QCD on a lattice provide increasingly accurate nu- 
merical values for these parameters, but little if any insight into their origin. The first approach relies 
on the intimate relation between QCD Green functions, in particular their Operator Product Expansion 
(OPE) beyond perturbation theory, and their hadronic counterparts. This relation follows from Cauchy's 
theorem in the complex energy plane, and is known as the QCD sum rule technique yy. In addition to 
producing numerical values for the QCD parameters, this method provides a detailed breakdown of the 
relative impact of the various dynamical contributions. For instance, the strong coupling at the scale of 
the T-lepton mass essentially follows from the relation between the experimentally measured t ratio, i?,-, 
and a contour integral involving the perturbative QCD (PQCD) expression of the V + A correlator. This 
is the cleanest, most transparent, and model independent determination of the strong coupling [2]-[3]. It 
also allows to gauge the impact of each individual term in PQCD, up to the currently known five-loop 
order. Similarly, in the case of the quark masses one considers a QCD correlation function which on 
the one hand involves the quark masses and other QCD parameters, and on the other hand it involves 
a measurable (hadronic) spectral function. Using Cauchy's theorem to relate both representations, the 
quark masses become a function of QCD parameters, e.g. the strong coupling, some vacuum condensates 



^Invited review paper to be published in Modern Physics Letters A. 



1 



reflecting confinement, etc., and measurable hadronic parameters. The virtue of this approach is that 

it provides a breakdown of each contribution to the final value of the quark masses. More importantly, 
it allows to tune the relative weight of each of these contributions by introducing suitable integration 
kernels. This last feature has been used recently to solve the historical problem of the systematic uncer- 
tainties affecting light quark mass determinations, to be discussed in this report. 

In the case of the light quark masses the ideal Green function is the pseudoscalar current correlator. This 
contains the square of the quark masses as an overall factor multiplying the PQCD expansion, and the 
leading power corrections in the OPE. Unfortunately, this correlator is not realistically accessible exper- 
imentally beyond the pseudoscalar meson pole. While the existence of at least two radial excitations of 
the pion and the kaon are known from hadronic interaction data, this information is hardly enough to 
reconstruct the full spectral functions. In spite of many attempts over the years to model them, there 
remains an unknown systematic uncertainty that has plagued light quark mass determinations from QCD 
sum rules. The use of the vector current correlator, for which there is plenty of experimental data from 
T decays and €^e~ annihilation, is not a realistic option for the light quarks as their masses enter as 
sub-leading terms in the OPE. The scalar correlator, involving the square of quark mass differences, at 
some stage offered some promise for determining the strange quark mass with reduced systematic un- 
certainties. This was due to the availability of data on if — tt phase shifts. Unfortunately, these data 
do not fully determine the hadronic spectral function. The latter can be reconstructed from phase shift 
data only after substantial theoretical manipulations, implying a large unknown systematic uncertainty. 
A breakthrough has been made recently by introducing an integration kernel in the contour integral in 
the complex energy plane. This allows to suppress substantially the unknown hadronic resonance con- 
tribution to the pseudoscalar current correlator. As it follows from Cauchy's theorem, this suppression 
implies that the quarks masses are determined essentially from the well known pseudoscalar meson pole 
and PQCD (well known up to five-loop level). In this way it has been possible to reduce the hadronic 
resonance contribution to the 1% level, allowing for an unprecedented accuracy of some 8 — 10 % in the 
values of the up-, down-, and strange-quark masses. Further improvement on this accuracy will be pos- 
sible with further reduction of the uncertainty in the strong coupling, now the main source of error. 
The determination of the charm- and bottom-quark masses has been free of systematic uncertainties due 
to the hadronic resonance sector, as there is plenty of experimental information in the vector channel 
from e^e~ annihilation into hadrons. One problem, though, is that the massive vector current correlator 
is not known in PQCD to the same level as the light pseudoscalar correlation function. Nevertheless, 
substantial theoretical progress has been made over the years leading to extremely accurate charm- and 
bottom-quark masses. The novel idea of introducing suitable integration kernels in Cauchy's contour inte- 
grals, as described above, has also been used recently as a way of improving accuracy in the heavy-quark 
sector. For instance, kernels can be used to suppress regions where the data is either not as accurate, or 
simply unavailable. This will also be reported here. 

The paper is organized as follows. First, determinations of quark-mass ratios from various hadronic data, 

as well as from chiral perturbation theory, will be reviewed in Section 2. These ratios are quite useful as 
consistency checks for results from QCD sum rules. Section 3 describes the OPE beyond perturbation 
theory, one of the two pillars of QCD sum rules. Section 4 discusses quark-hadron duality and finite 
energy sum rules. These sum rules weighted by suitable integration kernels will be analyzed in the light 
quark sector in Section 5. In particular it will be shown how this technique unveils the subjacent hadronic 
systematic uncertainty plaguing light quark mass determinations for the past thirty years. In Section 6 
recent progress on charm- and bottom-quark mass determinations will be reported. Comparison with 
lattice QCD results for all quark masses will also be made. Finally, Section 7 provides a very short 
summary of this report. 

As an important disclaimer, this paper is not a comprehensive review of past quark mass determinations 
from QCD sum rules. It is, rather, a report on recent progress on the subject. Given that past deter- 
minations of light quark masses were affected by unknown systematic uncertainties, essentially from the 
hadronic resonance sector, it makes little or no sense to review them once the main uncertainty has been 
exposed. Any agreement between values affected by this uncertainty and current results, free of it, would 
only be fortuitous. In fact, once this hadronic resonance uncertainty is removed, the values of all three 
light quark masses get reduced by some 15 - 20 %, a clear sign of a systematic uncertainty acting in 
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only one direction. Last but not least, light quark masses from QCD sum rules before 2006 employed 
correlators up to at most four-loop level in PQCD, together with superseded values of the strong coupling. 



2 Quark mass ratios 



Quark masses actually precede QCD by a number of years, albeit under the guise of current algebra quark 
masses, which clearly lacked today's detailed understanding of quark-mass renormalization. In fact, the 
study of global SU(3) x SU(3) chiral symmetry realized a la Nambu-Goldstone, and its breaking down to 
SU(2) X SU(2), followed by a breaking down to SU(2), and finally to U(l) was first done using the strong 
interaction Hamiltonian [l]-[7] 

H{x) = Ho{x) + eo uo{x) + £3 usix) + es usix) . (1) 

The term Ho{x) above is SU(3) x SU(3) invariant, the £0,3,3 are symmetry breaking parameters, and 
the scalar densities uo,3.8(a^) transform according to the 3 3 -1- 33 representation of SU(3) x SU(3). In 
modern language, eg is related to the strange quark mass rUs, and £3 to the difference between the down- 
and the up-quark masses md — , while the scalar densities are related to products of quark-anti-quark 
field operators. For instance, the ratio of SU(3) breaking to SU(2) breaking is given by 

rris- mud £3 

= : (,2) 

md-mu 2 £3 

where m^d = {mu + md)/2. In the pre-QCD era many relations for quark-mass ratios were obtained from 
hadron mass ratios, as well as from other hadronic information, e.g. 77 — >■ Stt, K13 decay, etc. [7]. To 
mention a pioneering determination of the ratio R above, from a solution to the 77 — >■ 37r puzzle proposed 
in [5] it followed [9, R^^ = 0.020 ± 0.002, in remarkable agreement with a later determination based on 
baryon mass splittings [TU] R^^ = 0.021 ± 0.003. With the advent of chiral perturbation theory (CHPT) 
[6]- [2]! |11|-|12|. certain quark mass ratios turned out to be renormalization scale independent to leading 
order, and could be expressed in terms of pseudoscalar meson mass ratios [7],llll-|13j, e.g. 



mu ^ M],+ - M|o + 2M^o - 
md Mlo - Ml^ + Ml^ 



0.56 (3) 



TO,, 



md Mlo - Ml+ + Ml 



20.2 , (4) 



1^0 - -r iyi^+ 

where the numerical results follow after some subtle corrections due to electromagnetic self energies [T2] . 
Beyond leading order in CHPT things become complicated. At next to leading order (NLO) the only 
parameter free relation is 

r,2 = TTil - ml^ ^ Ml - Ml Ml 
- ml-ml Ml„-Ml, M^ ' 

Other quark mass ratios at NLO and beyond depend on the renormalization scale, as well as on some 
CHPT low energy constants which need to be determined independently [II]-[T1]. After taking into 
account electromagnetic self energies, Eq.(5) gives [11] Q = 24.3, while a recent analysis of — 37r pT| . 
[I4] gives 

Q = 22.3 ±0.8. (6) 

The ratios R, Eq.(2), and Q, Eq.(5), together with the leading order ratios Eqs.(3)-(4), will prove useful 
for comparisons with QCD sum rule results. An additional useful quark mass ratio involving the ratios 
Eqs.(3)-(4) is 

TO, 2 m^/md , , 

= — ^ = = 28.1 ± 1.3 , (7) 

mud l + mu/md 

where the numerical value follows from the NLO CHPT relation [11], to be compared with the LO result 
from Eqs.(3)-(4), rg — 25.9, and a large Nc estimate [15] = 26.6 ± 1.6. 
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3 Operator product expansion beyond perturbation theory 



The OPE beyond perturbation theory in QCD, one of the two pillars of the sum rule technique, is an 
effective tool to introduce quark-gluon confinement dynamics. It is not a model, but rather a parametriza- 
tion of quark and gluon propagator corrections due to confinement, done in a rigorous renormalizable 
quantum field theory framework. Let us consider a typical object in QCD in the form of the two-point 
function, or current correlator 

niq^) ^ i f d^a;e'«" < 0| r( J(a;) J(0)) |0 >, (8) 



where the local current J{x) is built from the quark and gluon fields entering the QCD Lagrangian. 
Equivalently, this current can also be written in terms of hadronic fields with the same quantum num- 
bers. A relation between the two representations follows from Cauchy's theorem in the complex energy 
(squared) plane. This is often referred to as quark-hadron duality, the second pillar of the QCD sum 
rules method to be discussed in the next section. The QCD correlator, Eq.(8), contains a perturbative 
piece (PQCD), and a non perturbative one mostly reflecting quark-gluon confinement. The leading order 
in PQCD is shown in Fig.l. Since confinement has not been proven analytically in QCD, its effects can 
only be introduced effectively, e.g. by parameterizing quark and gluon propagator corrections in terms 
of vacuum condensates. This is done as follows. In the case of the quark propagator 

1 1 

Sf{p) = -7 =^ -7 , 2\ ' (^) 

the propagator correction contains the information on confinement, a purely non-perturbative 

effect. One expects this correction to peak at and near the quark mass-shell, e.g. for p ~ in the case 
of light quarks. Effectively, this can be viewed as in Fig. 2, where the (infrared) quarks in the loop have 
zero momentum and interact strongly with the physical QCD vacuum. This effect is then parameterized 
in terms of the quark condensate (0|g(0)q(0)|0). Similarly, in the case of the gluon propagator 



Fi gure ll Leading order PQCD correlator. All values of the four- momentum of the quarks in the loop are allowed. The 
blue wiggly Une represents the current of momentum q (—q^ » 0). 



the propagator correction A(fc^) will peak at fc ~ 0, and the effect of confinement in this case can be 
parameterized by the gluon condensate (0|as G^'^ ■ G^u\^) (see Fig. 3). In addition to the quark and 
the gluon condensate there is a plethora of higher order condensates entering the OPE of the current 
correlator at short distances, i.e. 

Coi C2^^+2(g^A^') (0|O2jv+2(a*')|0) , (11) 

where /x^ is the renormalization scale, and where the Wilson coefficients in this expansion, C2N +2{<f' , M^): 
depend on the Lorentz indices and quantum numbers of J{x) and of the local gauge invariant operators 
On built from the quark and gluon fields. These operators are ordered by increasing dimensionality and 
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the Wilson coefficients, calculable in PQCD, fall off by corresponding powers of — g^. In other words, this 
OPE achieves a factorization of short distance effects encapsulated in the Wilson coefficients, and long 
distance dynamics present in the vacuum condensates. Since there are no gauge invariant operators of 
dimension d = 2 involving the quark and gluon fields in QCD, it is normally assumed that the OPE starts 
at dimension d = 4. This is supported by results from QCD sum rule analyses of r-lepton decay data, 
which show no evidence of d = 2 operators [l6]-[T7]. The unit operator / in Eq.(ll) has dimension d = 




Figure 2: Quark propagator modification due to (infrared) quarks interacting with the physical QCD vacuum, and 
involving the quark condensate. Large momentum flows through the bottom propagator. 



and CqI stands for the purely perturbative contribution. The Wilson coefficients as well as the vacuum 
condensates depend on the renormalization scale. For light quarks, and for the leading d = 4 terms 
in Eq.(ll), the /x^ dependence of the quark mass cancels the corresponding dependence of the quark 
condensate, so that this contribution is a renormalization group (RG) invariant. Similarly, the gluon 
condensate is also a RG invariant, hence once determined in some channel these condensates can be used 
throughout. The numerical values of the vacuum condensates cannot be calculated analytically from 
first principles as this would be tantamount to solving QCD exactly. One exception is that of the quark 
condensate which enters in the Gell-Mann-Oakes-Renner relation, a QCD low energy theorem following 
from the global chiral symmetry of the QCD Lagrangian [TS]. Otherwise, it is possible to extract values 
for the leading vacuum condensates using QCD sum rules together with experimental data, e.g. e'^e~ 
annihilation into hadrons, and hadronic decays of the r-lepton. Alternatively, as lattice QCD improves 
in accuracy it should become a valuable source of information on these condensates. 




Figure 3: Gluon propagator modification due to (infrared) gluons interacting with the physical QCD vacuum, and involving 
the gluon condensate. Large momentum flows through the quark propagators. 



4 Quark- hadron duality and finite energy QCD sum rules 

Turning to the hadronic sector, bound states and resonances appear in the complex energy (squared) 
plane (s-plane) as poles on the real axis, and singularities in the second Riemann sheet, respectively. All 
these singularities lead to a discontinuity across the positive real axis. Choosing an integration contour as 
shown in Fig. 4, and given that there are no other singularities in the complex s-plane, Cauchy's theorem 
leads to the finite energy sum rule (FESR) 

f ° ds- fis) Imnis)\HAD ^ i ds J{s)n{s)\QCD , (12) 

Jsth ^TTZ Jc(\so\) 
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Im(s) 




Figure 4: Integration contour in the complex s-plane. The discontinuity across the real axis brings in the hadronic spectral 
function, while integration around the circle involves the QCD correlator. The radius of the circle is sqi the onset of QCD. 



where /(s) is an arbitrary (analytic) function, sth is the hadronic threshold, and the finite radius of the 
circle, sq, is large enough for QCD and the OPE to be used on the circle. Physical observables determined 
from FESR should be independent of sq. In practice, though, this is not exact, and there is usually a 
region of stability where observables are fairly independent of so, typically somewhere inside the range 
So — 1 — 4 GeV^. Equation (12) is the mathematical statement of what is usually referred to as quark- 
hadron duality. Since QCD is not valid in the time-like region (s > 0), in principle there is a possibility 
of problems on the circle near the real axis (duality violations), to be discussed shortly (this issue was 
identified very early in 19] long before the present formulation of QCD sum rules). The right hand side 
of this FESR involves the QCD correlator which is expressed in terms of the OPE as in Eq.(ll). The left 
hand side involves the hadronic spectral function which is written as 

I'mI{{s)\HAD = ImI{{s)\poLE + ImI{{s)\REsO{sQ- s) + ImIV{s)\pQCD 0{s - sq) , (13) 

where the ground state pole is followed by the resonances which merge smoothly into the hadronic 
continuum above some threshold so. This continuum is expected to be well represented by PQCD if so 
is large enough. Hence, if one were to consider an integration contour in Eq.(12) extending to infinity, 
the cancellation between the hadronic continuum on the left hand side and the PQCD contribution on 
the right hand side, would render the sum rule a FESR. The integration in the complex s-plane of the 
QCD correlator is usually carried out in two different ways. Fixed Order Perturbation Theory (FOPT) 
and Contour Improved Perturbation Theory (CIPT). The first method treats running quark masses and 
the strong coupling as fixed at a given value of so- After integrating all logarithmic terms (ln(— s//i^)) 
the RG improvement is achieved by setting the renormalization scale to /i^ = — so. In CIPT the RG 
improvement is performed before integration, thus eliminating logarithmic terms, and the running quark 
masses and strong coupling are integrated around the circle. This requires solving numerically the RGE 
for the quark masses and the coupling at each point on the circle. The FESR Eq.(12) with /(s) = 1 and 
in FOPT can be written as 



If the hadronic spectral function is known in some channel from experiment, e.g. from r-decay into 
hadrons, then ImIi{s)\HAD = I''nIl{s)\DATA, and Eq.(14) can be used to determine the values of the 




2Ar+2(so) , 



(14) 



where the dimensionless PQCD moments M2n+2{so) are given by 




(15) 
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Figure 5: Results for f„ from the standard FESR in the axial-vector channel, Eq. (14) with N=0, with no dimension 
d = 2 term, and using CIPT, with Aqcd = 365 MeV (as(MT-) = 0.335). The straight line is the experimental value of 
and the points are the integrated data with the experimental errors. 



vacuum condensates. Subsequently, Eq.(14) can be used in a different channel for a different application. 
It is important to mention that the correlator n((7^) is generally not a physical observable. However, this 
has no effect in FOPT as the unphysical quantities (polynomials) in the correlator do not contribute to 
the integrals. In the case of CIPT, though, this requires modified sum rules involving as many derivatives 
of the correlator as necessary to render it physical. 

Next, let us consider an application where the integration kernel /(s) in Eq.(12) is of great importance 
[20| . In the axial- vector channel, the FESR Eq.(14) with /(s) — 1 and = can be confronted with 
data from r-decay. The hadronic spectral function is then written as the sum of the pion pole and 
the resonance data known up to the kinematical end point sq = M^. The moment M2{sq) is known 
up to five-loop order in PQCD, so that the FESR can be used to confront the resonance data plus 
PQCD with e.g. /^. As seen from Fig. 5 the agreement is rather poor, except possibly near the end 
point. At first sight, this may be interpreted as a signal for quark-hadron duality violations near the 
real axis, even at this high enough energy. In fact, it has been known for quite some time that e.g. 
the Weinberg (chiral) sum rules are not saturated by the t decay data unless one introduces pinched 
integration kernels, e.g. /(s) = [1 — (s/so)]^^+^^ Unfortunately, the r-lepton is not massive 

enough to probe higher energy regions. In spite of this it is still possible to explore a wider energy range 
by introducing as integration kernel a polynomial /(s) = P(s,so,si) tuned to eliminate the (unknown) 
hadronic contribution to the integral between si and sq> si, where si is at or near the end point of the 
data. It has been shown [20] that the optimal degree of P{s) is the simplest, i.e. the linear function 



P(s, so,si) = 1 - 



2s 



So + si 



so that 



constant x / P{s, sq: si)ds = . 

J si 



(16) 



(17) 



In this case the complete FESR becomes a linear combination of a dimension-two and a dimension-four 
FESR, which from Eqs.(14) and (16) is given by 



2f = - 



[ ' dsP{s)-Im\l{s)\DATA + 



M2(so) - ^^M4(so) 



So + Si 



+ 



1 

4^ 



^2(02 



So + Si 



A(so) 



(18) 
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Figure 6: Results for from the FESR in the axial-vector channel, Eq. (18), with 02(62) = 0, 04(64) = 
using CIPT with Aqcd = 365 MeV (asiMr) = 0.335). The straight line is the experimental value of /„ 
with errors are the integrated data up to ~ Mt- 



: 0.05 GeV*, and 
, and the points 



where the pion pole has been separated from the data, and the chiral hmit is understood. The term 
A(so) is the error being made by assuming that the data is constant in the interval si — sq- It is possible 
to estimate this error which turns out to be two to three orders of magnitude smaller than 2/^ on the left 
hand side of Eq.(18) [20]. As can be seen from Fig. 6 the FESR Eq.(18) shows an excellent consistency 
between QCD and the r data in the axial- vector channel in a remarkably wide region sq — 4 — 10 GeV^. 
A similar consistency is also found in the vector channel, where QCD is now confronted with zero (there 
is no pole in this channel). This result shows either no evidence for quark-hadron duality violations in 
these channels, or if they are present it indicates that they are suppressed by the integration kernel (some 
model dependent analyses claim the existence of duality violations ^24| ) . 



5 Light quark masses 



Traditionally, the light quark masses have been determined using the correlator, Eq.(8), involving the 
pseudoscalar currents J{x) = dij,A^{x)\j — [mi(/i) +mj{^)] : qj{x)ij5qi{x) :, where A^{x) is the axial 
vector current of flavours i and j, Wiidi) the quark mass in the MS scheme, /i the renormalization scale 
and qi{x) are the quark fields. An issue of major concern in the past was the presence of logarithmic 
quark-mass singularities in these correlators. This problem has been satisfactorily resolved some time 
ago in [IS]- [55]. These correlators are now known to five-loop order in PQCD [57], and free of logarithmic 
quark mass singularities. The Wilson coefficients of the leading power corrections, i.e. the gluon and the 
quark condensates, are also known up to two-loop level [25]-[26]. Higher dimensional condensates, as well 
as quark mass corrections of order 0{mf) (with respect to the one- loop term) and higher turn out to be 
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negligible. From Cauchy's theorem, Eq.(12), the FESR to determine the quark masses can be written as 



27^* Jci\so\) 



rso I 

/ ds- Im i^5{s)\RES A5(s) , (19) 



where A5(s) is an (analytic) integration kernel to be introduced shortly, the first term on the right 
hand side is the pseudoscalar meson pole contribution (P — tt.K), sth is the hadronic threshold, and 
Im'iJj5{s)\RES is the hadronic resonance spectral function. The radius of integration so is assumed to be 
large enough for QCD to be valid on the circle. For later convenience this FESR can be rewritten as 

<55(so)|Qcr> = S^Ipole + S5{sq)\b^eS , (20) 

where the meaning of each term is self evident. Historically, the problem with the pseudoscalar correlator 
has been the lack of direct experimental information on the hadronic resonance spectral functions. Two 
radial excitations of the pion and of the kaon, with known masses and widths, have been observed in 
hadronic interactions |28j . However, this information is hardly enough to reconstruct the full spectral 
function. In fact, inelasticity, non-resonant background and resonance interference are impossible to guess, 
leaving no choice but to model these functions. This introduces an unknown systematic uncertainty which 
has been present in all past QCD sum rule determinations of the light quark masses. Since the FESR 
Eq.(19) is valid for any analytic A5(s) one can choose this kernel in such a way as to suppress (55(so)|i?£;s 
as much as possible. An example of such a function is the second degree polynomial |^-|32) 

A5(s)|fl.Es = 1 - oo s - ai , (21) 

where oq and ai are constants fixed by the requirement A5(Aff) = A5(Af|) = 0, where Mi. 2 are the 
masses of the first two radial excitations of the pion or kaon. This simple kernel suppresses enormously 
the resonance contribution, which becomes only a couple of a percent of the pole contribution, and well 
below the current uncertainty due to the strong coupling. This welcome feature is essentially independent 
of the model chosen to parametrize the resonances. A practical parametrization consists of two Breit- 
Wigner forms normalized at threshold according to chiral perturbation theory, as first proposed in [33j 
for the pionic channel, and in |34] for the kaonic channel. Detailed results for S5{so)\qcd, to five-loop 
order in PQCD and up to dimension d = 4 in the OPE, after integrating in FOPT may be found in [50] . 
In the case of CIPT the FESR must be written in terms of the second derivative of the current correlator. 
This is in order to eliminate the unphysical first degree polynomial present in ipbis), which unlike the 
case of FOPT would otherwise contribute to the FESR which then becomes 

ds ^;'<3^^(s)[^^(s)-F(so)]=2/|A/|A5(Af|) 

27n Jcdsnl) 



So 



ds lTmP5{s)\REs ^5is) , (22) 

Sth 



where 

^'0 „ *o 



and 



F{s) = ^.s[so~aof-a,^\+--aoj^a,— , (23) 



„2 3 „4 

0>n -^n 



F(.o) = -f + ao|+ai^. (24) 

The RG improvement is used before integration, so that all logarithmic terms vanish. The running 
coupling as well as the running quark masses are no longer frozen as in FOPT, but must be integrated. 
This can be done by solving numerically the respective RG equations at each point on the integration 
circle in the complex s-plane. Detailed expressions are given in |30)-|31j. 
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Fi gurc 7l The strange quark mass ms{2 GeV) in the MS scheme taking into account only the kaon pole with A5{s) = 1 
(curve(a)), and the two Breit-Wigner resonance spectral function with a threshold constraint from CHPT |34) . with A5(s) = 
1 (curve(b)), and A5{s) as in Eq.(21) (curve(c)). A systematic uncertainty of some 20% due to the resonance sector is 
dramatically unveiled. 



The parameters of the integration kernel, Eq.(21), are oq — 0.897 GeV~^, and ai — —0.1806 GeV"'' 
for the pionic channel, and Gq = 0.768 GeV"^, and oi — —0.140 GeV"'* for the kaonic channel. These 
values correspond to the radial excitations 7r(1300), 7r(1800), K{1A60) and A'(1830). The pion and kaon 
decay constants are ^ U = 92.21 ± 0.14 MeV, and /k = (1.22 ± 0.01)/^. In the QCD sector it is best 
to use the value of the strong coupling determined at the scale of the r-mass, as this is close to the scale 
in current use for the light quark masses, i.e. /i = 2 GeV. The extraction of as{Mr) from the Rr ratio 
involves an integral with a natural kinematical integration kernel that eliminates the contribution of the 
d = 4 term in the OPE. This welcome feature improves the accuracy of the determination, and it makes 
little sense to introduce additional spurious integration kernels which would artificially recover this d — A 
contribution. The different values obtained from r decay using CIPT are all in agreement with each 
other, i.e. a,{Mr) = 0.338 ± 0.012 gj, a^(M^) = 0.341 ± 0.008 a^iMr) = 0.344 ± 0.009 [371, and 
as{Mr) = 0.332 ± 0.016 [38]. These determinations are model independent and extremely transparent, 
with as obtained essentially by confronting PQCD with the single experimental number Rt-. The method 
of FOPT is known to give rise to a pathological non-convergent perturbative series in this application 
[2] , so it will not be considered here. The d — 4 gluon condensate has been extracted from r decays [17] , 
but one can conservatively consider the wide range< a^G^ >= 0.01 — 0.12 GeV*. The impact of the 
light quark condensate is at the level of 1% in the quark masses. A ± 30% uncertainty in the resonance 
contribution (55(so)|-R-ES in Eq.(20) translates into a safe 1% change in the quark masses. Finally, it has 
been assumed that the unknown six-loop PQCD contribution is equal to the five-loop result, an extreme 
but very conservative estimate of higher orders in PQCD. 

Beginning with the strange quark mass. Fig. 1 shows the results for ms(2 GeV)\jjg- with no integration 
kernel, A5(s) = 1, and taking into account only the kaon pole, curve (a), and the kaon pole plus a two 
Breit-Wigner resonance model with a threshold constraint from CHPT [34], curve (b) (a misprint in the 
formula for the spectral function in |34j has been corrected in [35]). These curves are for the central value 
of as{Mr) whose uncertainties will be considered afterwards. The latter result is reasonably stable in the 
wide region sq = 2 — 4 GeV^, so that it could lead us to conclude that ms(2 GeV)|jg^ ~ 100 — 120 MeV, 
albeit with a yet unknown systematic uncertainty arising from the resonance sector. Introducing the 
kernel, Eq.(21), leads to curve (c) and to a dramatic unveiling of this systematic uncertainty. In fact, the 
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Source 


TUu 




ms 


mud 


mu/md 


ms/mud 


R 





QCDSR [31] 


2.6 ±0.3 


5.6 ±0.4 




4.1 ±0.3 


0.46 ±0.06 








QCDSR [30] 






102 ±8 






24.9 ± 2.7 


33 ±6 


21 ± 3 


FLAG [S] 


2.2 ± 0.3 


4.6 ±0.6 


95 ±10 


3.4 ±0.4 


0.47 ± 0.04 


27.8 ± 1.0 


37.2 ±4.1 


23.1 ± 1.5 



Table 1: The running quark masses in the MS scheme at a scale ^ = 2 GeV in units of MeV from QCD sum rules 
(QCDSR) (first two rows), and from the FLAG lattice QCD analysis 1121 . The ratios Wiu/fnd Eind m,s/rn^d a-re an input 
in the QCDSR (see text). The ratios R and Q are defined on the l.h.s. of Eqs.(2) and (3). The QCDSR results for the up- 
and down- quark masses and ratios differ slightly from |31 | -|32 | due to the input value oimu/rrid used here. 

real value of the quark mass is ffis (2 GeV) = 102 ± 8 MeV, or some 20% below the former result (this 
error now includes the uncertainty in as). In addition, and as a bonus the systematic uncertainty- free 
result is remarkably stable in the unusually wide region sq — 1 ~ 4 GeV^ (typical stability regions are 
only half as wide) . 

It must be recalled that the pseudoscalar correlator involves the overal factor (m^ + niud)^ ■ Hence, in 
order to determine nis an input value for the ratio ms/niud is needed in the result from the sum rule, 
which is 

, ^ 105.5 ±8.2 MeV , , 

m.(2 GeV)|jg= . (25) 

1 ± "Itid/'Tls 

Using the wide range ms/rriud = 24 — 29 leads to Wis (2 GeV)\jj-^ = 102 ± 8 MeV. In this case the impact 
of the uncertainty in the quark mass ratio is small. However, in the case of the up- and down-quark 
masses the corresponding ratio mu/md plays a more important role in the result from the sum rule, 
which is 

8.2 ±0.6 MeV 

1 ± mu/md 

The input used in Eq.(26) for the ratio mu/md = 0.47 ± 0.04 is from the overall lattice QCD analysis 
of the FLAG group [T^j. Once is determined from Eq.(26), m„ follows. Using these results for 
the individual masses one obtains the ratios mu/md and ms/mud shown in Table 1. The quark masses 
mu(2 GeV)\jjg and md{'2 GeV)|p^ also exhibit a remarkably wide stability region sq ~ 1 — 4 GeV^ 
[21]- [35]. In Tabic 1 one finds a summary of the results for the light quark masses, and the ratios R and 
Q defined on the left hand side of Eqs. (2) and (5), together with the results of the Flag group [T^ . 
The values of the up- and down-quark masses and their ratios in Table 1 are slightly different from those 
in ^31]-[32 due to the input value for the ratio mu/md (in [31j-|32| the value mu/md = 0.55 was used). 
The various sources of errors in the quark masses discussed earlier combine into the final values given in 
Table 1. Having all but eliminated the systematic uncertainty from the hadronic resonance sector, the 
main source of error is now due to the strong coupling. Improved accuracy in the determination of as 
would then allow for a reduction of the uncertainties in the light quark masses. 



md{2 GeY)\j^ = . (26) 



6 Heavy quark masses 

Determinations of the charm- and bottom-quark masses are not affected by a lack of data, as there is 
plenty of experimental information from e^e~ annihilation into hadrons |28) . except for a gap in the 
region 25 GeV^ ^ ^ < 50 GeV^ . On the theoretical side there has been very good progress on PQCD up 
to four- loop level 39^ -[5F. The leading power correction in the OPE is due to the gluon condensate with 
its Wilson coefficient known at the two- loop level [56j- The correlator, Eq.(8), involves the vector current 
J{x) = Vi^{x) = Q{x)^fj,Q{x), where Q{x) is the charm- or bottom-quark field. The experimental data 
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is in the form of the i?Q-ratio for charm (bottom) production, which determines the hadronic spectral 
function. Modern determinations of the heavy-quark masses have been based on inverse moment (Hilbert- 
type) QCD sum rules, e.g. Eq.(12) with f{s) — 1/s". These sum rules require QCD knowledge of the 
vector correlator in the low energy region, around the open charm (bottom) threshold, as well as in the 
high energy region. A recent update [57] of earlier determinations [39]- [41], [43]- [45] reports a charm- 
quark mass in the MS scheme accurate to 1%, and half this uncertainty for the bottom-quark mass. 
However, the analysis of |55| claims an error a factor two larger for the charm-quark mass. It appears 
that the discrepancy arises from the treatment of PQCD. In fact, in [5S] two different renormalization 
scales were used, one for the strong coupling and another one for the quark mass. This unconventional 
choice results in a much larger error in the charm-quark mass obtained from inverse (Hilbert) moment 
QCD sum rules. It does not affect, though, sum rules involving positive powers of s. In any case, the 
philosophy in current use is to choose the result from the method leading to the smallest uncertainty. 
Beginning with the charm-quark mass, an alternative procedure was proposed some years ago based only 
on the high energy expansion of the heavy-quark vector correlator [5S]-[Sn]. This method was followed 
recently [60], but with updated PQCD information and the inclusion of integration kernels in the FESR, 
Eq.(12), tuned to enhance/suppress contributions from data in certain regions. The first such kernel is 
the so-called pinched kernel [21)-[^ 

f{s) = l--, (27) 
So 

which is supposed to suppress potential duality violations close to the real s-axis in the complex s-plane. 
In connection with the charm-quark mass application, this kernel enhances the contribution from the 
first two narrow resonances, J/tp and tp{2S), and reduces the weight of the broad resonance region, 
particularly near the onset of the continuum. The latter feature is better achieved with the alternative 
kernel [HT] 



/(^) = l-(fj , (28) 

which produces an obvious larger enhancement of the narrow resonances, and a larger quenching of the 
broad resonance region. This kernel has been used together with both the high and the low energy 
expansion of the vector correlator in [61]. Since there are no data in the region 25 GeV^ ^ s < 50 GeV^, 
while the data for s > 50 GeV^ agrees with PQCD, it is useful to introduce a kernel that will allow a 
suppression in the former region, e.g. 

/(s) =P„[x(s)] , (29) 

where 

.(.) = '^'(^° + -^^) , (30) 
So - si 

with So > si, and Vn{x) are the standard Legendre polynomials, i.e. Viix) = x, V2{x) — {5x^ — 3a::)/2, 
etc., which satisfy the constraint 

/ s'' Vn[x{s)] ds = , (31) 

J Si 

where si ~ 24 GeV^, and sq varies in the region where there is no data. These Legendre- type kernels 
provide extra weight to the well known resonance region on account of their rapid growth for s < si. The 
resulting charm-quark masses are essentially insensitive to the order of these polynomials, with n — 3 — 6 
giving answers differing at the 0.1% level. 

The results for TOc(3GeV) in the MS scheme using various integration kernels are listed in Table 2, 
together with the various uncertainties. The merits of each kernel may be judged by its ability to minimize 
these uncertainties, in particular those that might be most affected by systematic errors, such as e.g. the 
experimental data. The kernel /(s) = 1 — (so/s)^ appears to be optimal as it produces the smallest 
uncertainty due to the data, and is very stable against changes in so- Some recent determinations of 
the charm-quark mass are based on Hilbert moments with no Sg-dependent kernel, corresponding to the 
third row in Table 2. While there is no explicit SQ-dependence in Hilbert moments (the integrals extend 
to infinity), there is definitely a residual dependence when choosing the threshold for the onset of PQCD. 
From a FESR perspective, the major drawback of the kernel 1/s^ is clearly the poor stability against 
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Uncertainties (in MeV) 



Kernel 


rnci3 GeV) 


EXP 


Aas 




NP 


so 


Total 


1 - s/so 


983 


9 


1 




1 


16 


25 


nix{s)] 


1007 


22 


1 


8 


2 


< 1 


23 




995 


9 


3 


1 


1 


14 


17 


1 - {so/s)^ 


987 


7 


4 


1 


1 


4 


9 



Table 2: Results for mc(3GeV) in the MS scheme and in MeV from different integration kernels in the FESR. Legendre 
polynomial kernels for n = 3 — 6 give basically the same result as for n = 5. The various uncertainties are due to the data 
(EXP), the value of as (Aa^), changes of ±35% in the renormalization scale around fi = 3 GeV (A/i), the value of the 
gluon condensate (NP), and due to variations of sg (•^o) in a very wide range (see text). The first four errors are added in 
quadrature, but the last one is conservatively added linearly. The first two rows are from 1601 and the last two from I61| . 
The total error in the first row includes an estimate of the error due to PQCD truncation at the four-loop level (2 MeV). 
CIPT was used throughout. 



changes in sq. In this regard, the Legendre polynomial kernels would be optimal, except that they have 
the largest uncertainty due to the data. An important remark is in order concerning the uncertainty 
due to changes in sq. From current data it is not totally clear where does PQCD start. This problem 
not only affects FESR, with their explicitly obvious so-dependence, but also Hilbert moments with an 
implicit SQ-dependence, as there is no data all the way up to infinity. The range of variation in sq and 
its contribution to the uncertainty in the quark mass, as appearing in Table 2, is as follows. For rows 1, 
3 and 4 sq — 15 — 23 GeV^, and for the Legendre polynomial kernels, row 2, it is sq — 100 — 200 GeV^. 
Two of the most recent results for mc(3 GeV) [57], [55], together with the weighted FESR value [61] (last 
row in Table 2) are 

r 986 ± 13 MeV [57] 

mc(3GeV) = < 998 ± 29 MeV [55] (32) 
I 987 ± 9 MeV [M] , 

in very good agreement with each other, except for the errors. The small uncertainty from [61] is due in 
part to improved quenching of the data in the broad resonance region, but mostly due to a strong reduction 
in the sensitivity to sq, i.e. the onset of PQCD. For comparison, a recent lattice QCD determination 
gives [62] 

TOc(3 GeV) = 986 ± 6 MeV . (33) 

Turning to the bottom-quark mass, the most recent determination from Hilbert moment QCD sum rules 
[57] is 

r?T6(10GeV) = 3610 ± 16 MeV . (34) 

The error above is due to uncertainties in the data, the strong coupling and the renormalization scale. It 
does not include, though, an uncertainty due to the onset of PQCD. There is an appreciable difference 
if PQCD were to start at the end point of the (BABAR) data {^/sq ~ 11.2 GeV), or rather at higher 
energies (-yio — 13.0 GeV). If one were to adopt the procedure of [6^, i.e. add in quadrature the first 
three errors, but only linearly the one due to so, then using the information in |57j the result would 
change to TOf,(10 GeV) = 3600 ± 47 MeV. Presumably, the error would be larger if a broader region 
in So were to be considered. One of the virtues of weighted sum rules is to reduce this uncertainty, as 
may appreciated from Table 2 for the case of the charm-quark mass. For the bottom-quark mass work 
on weighted FESR is in progress [63]. In any case, it should be kept in mind that by choosing different 
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kernels it is possible to generate a large number of predictions for the quark masses with a variety of 
different errors. According to current philosophy one chooses the determination having the smallest error. 



7 Conclusions 

After a short review of quark mass ratios the method of QCD sum rules was discussed, in connection with 
determinations of individual values of the quark masses. The historical (unknown) hadronic systematic 
uncertainty affecting light quark mass determinations was highlighted. Details of the recent breakthrough 
in eliminating this uncertainty were provided. Future improvement in accuracy is now possible, and 
depends essentially on more accurate determinations of the strong coupling, now the main source of 
error. The new values of the light quark masses and their ratios, free from this systematic uncertainty, 
agree well with lattice QCD results. In the heavy quark sector recent more accurate determinations of 
the charm- and bottom-quark masses were reported. While these values are all in agreement, there is 
some disagreement on the size of the errors. The use of suitable multipurpose integration kernels in 
FESR allows to tune the weight of the various contributions to the quark masses. This in turn allows 
to minimize the error due to the data, as well as to the uncertainty in the onset of PQCD. The latter 
uncertainty impacts FESR as well as Hilbert moment sum rules, as there is no data all the way up to 
infinity. If no kernel, other than inverse powers of s is used then this uncertainty would be much larger 
than normally reported, as may be appreciated from Table 2. However, according to current philosophy 
one chooses the determination having the smallest error. Marginal improvement of the current total error 
in this framework should be possible with improved accuracy in the data and in the strong coupling. 
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